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ABSTRACT 


An  error  probability  analysis  is  performed  for  a  conventional  noncoherent  M-ary 
orthogonal  frequency-shift  keying  (MFSK)  receiver  employing  fast  frequency-hopped 
(FFH)  spread  spectrum  waveforms  transmitted  over  a  frequency-nonselective,  slowly 
fading  Ricean  channel  with  partial-band  noise  interference  .  Each  diversity  reception  is 
assumed  to  fade  independently.  The  partial-band  interference  is  modeled  as  a  Gaussian 
process.  The  effects  of  wideband  thermal  noise  are  also  included.  The  energy  per  hop  is 
held  constant;  thus,  as  diversity  increases,  energy  per  symbol  increases.  Previous 
analyses  considered  only  constant  energy  per  symbol  systems,  however,  practical  military 
systems  are  likely  to  employ  fixed  hop  rates.  There  is  some  performance  enhancement  to 
be  obtained  from  implementing  diversity  in  a  conventional  FFH/MFSK  system  with  fixed 
hop  rates,  but  partial  band  interference  still  results  in  significant  degradation. 

Additionally,  the  performance  of  this  FFH  receiver  is  investigated  over  the  same 
channel  in  the  presence  of  partial-band  tone  jamming  without  diversity  for  the  case  of 
binary  frequency-shift  keying  (BFSK)  when  both  the  signal  and  the  jammer  can  fade 
independently.  Performance  when  only  a  single  jamming  tone  per  hop  slot  is  allowed  is 
compared  to  that  obtained  when  two  jamming  tones  per  hop  slot  arc  possible.  When  the 
jamming  signal  experiences  Rayleigh  fading  there  is  very  little  degradation  of  the 
jammer's  effectiveness  as  compared  to  when  the  jamming  signal  is  not  affected  by  fading. 
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I.  INTRODUCTION 


A  fast  frequency-hopping  (FFH)  communication  system  is  a  subset  of  spread 
spectrum  communications  that  utilizes  a  bandwidth  greatly  exceeding  that  required  for 
the  information  signal  alone.  Frequency-hopping  spread  spectrum  is  fundamentally 
different  from  direct  sequence  (DS)  spread  spectrum  in  the  technique  of  signal  generation 
and  recovery.  This  thesis  focuses  on  frequency-hopping  spread  spectrum  systems 
because  of  their  practical  military  importance.  As  spread  spectrum  systems  grow  more 
popular  and  occupy  wider  communications  bandwidth,  the  likelihood  of  both  hostile  and 
non-hostile  sources  of  narrowband  interference  also  grows.  It  is  important  for  mission 
planners  to  have  a  reliable  estimate  of  the  degradation  that  their  communication  systems 
will  suffer  under  partial-band  noise  interference  as  well  as  tonal  interference.  Fast 
frequency-hopping  spread  spectrum  techniques  have  evolved  to  counter  the  threat  of 
intentional  jamming  [Ref.  1  ].  This  thesis  presents  an  error  probability  analysis  of  the 
conventional  fast  frequency-hopped  orthogonal  M-ary  frequency-shift  keying 
(FFH/MFSK)  receiver  with  noncoherent  detection  for  communications  over  channels 
with  Ricean  fading  of  the  signal  and  partial-band  noise  interference.  The  influence  of 
partial-band  tonal  interference  on  the  binary  orthogonal  frequency-shift  keying 
(FFH/BFSK)  receiver  without  diversity  where  both  the  signal  and  jammer  may  fade 
independently  is  also  considered. 
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A.  FFH/MFSK 


Frequency-shift  keying  is  popular  as  a  signaling  scheme  because  it  allows  for 
noncoherent  reception  of  the  signal.  A  typical  MFSK  signal  set  can  be  expressed  as 

si(0  =  <j2Accos(2-nf0t  +  Q\)  (1) 

s2(t)=  j2Accos[2n(f0  +  Af)t  +  Q2]  (2) 

and  so  on  to 

sa/(0=  j2Accos[2n(f0  +  (M-  1)A f)t  +  QM]  (3) 

where/,  is  the  lowest  signal  tone,  A/is  the  tone  spacing,  and  0m  is  the  phase  associated 
with  each  tone.  This  set  is  then  modulated  by  a  carrier  that  varies  pseudorandomly  in 

frequency  for  transmission  as  a  frequency-hopped  signal. 

All  communication  signals  suffer  from  interference  and  noise.  One  method  to 
overcome  this  degradation  is  to  transmit  the  signal  more  than  once  thus  providing  a  form 
of  diversity.  Fast  frequency-hopping  employs  this  redundancy  as  well  as  deliberately 
spreading  the  bandwidth.  Frequency-hopping  systems  in  which  several  symbols  are 
transmitted  per  hop  are  considered  slow  frequency-hopping,  while  those  transmitting 
several  hops  per  symbol  are  considered  fast  frequency-hopping.  Our  FFH/MFSK 
transmitter  performs  L  hops  per  data  symbol,  which  results  in  a  diversity  of  L.  At  the 
receiver  the  dehopped  signals  are  recovered  noncoherent!)  by  two  correlators  in  phase 
quadrature  with  the  dehopped  signal  waveform.  The  correlator  outputs  are  sampled  even 
Th  seconds,  where  Th  is  the  hop  period.  The  sampled  output  of  each  correlator  pair  is 
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squared  and  then  these  outputs  are  summed  L  times  to  obtain  the  decision  statistic  for 
each  branch  of  the  M-ary  detector.  The  largest  signal  detected  is  selected  as  the 
transmitted  symbol.  A  typical  receiver  structure  is  shown  in  Figure  1 . 

Proper  reception  and  demodulation  of  the  spread  signal  depend  on  the  recovery  of 
several  pieces  of  timing  information.  Both  the  sender  and  receiver  need  the  same 
pseudorandom  sequence  operating  synchronously.  Also  required  are  the  symbol  period 
and  the  hop  period.  In  practice  these  are  estimated  from  the  received  waveforms.  This 
thesis  assumes  that  this  information  is  recovered  without  error. 

B.  MULTI-PATH  EFFECTS 

The  losses  experienced  by  the  signal  during  propagation  is  worthy  of  an  entire  study 
in  itself.  However,  it  is  useful  to  make  some  general  observations  on  the  composition  of 
the  received  signal.  It  is  possible  and  even  likely  that  the  received  signal  arrives  at  the 
receiver  after  transiting  a  variety  of  different  paths.  Signals  traveling  a  longer  distance 
arrive  delayed  relative  to  the  direct  path  signal.  This  leads  to  multi-path  effects.  The 
magnitude  of  the  multi-path  effects  depends  on  the  magnitude  of  the  delayed  signal 
strength  versus  the  direct  path  signal  strength.  It  is  common  to  consider  the  sum  of  all 
delayed  signals  as  the  diffuse  component  of  the  received  signal.  [Ref.  2] 
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The  consequence  of  multi-path  reception  is  to  cause  the  signal  to  fade  in  a  time 
varying  fashion.  We  can  broadly  characterize  the  channel  conditions  by  examining  the 
magnitude  of  the  direct  signal  power  to  the  diffuse  signal  power.  The  Ricean  channel  is 
the  general  case.  Channels  that  have  nearly  all  the  received  signal  energy  in  the  direct 
component,  i.e.,  direct-to-diffuse  ratio  greater  than  100,  have  essentially  no  fading.  In  the 
limit,  an  infinite  direct-to-diffuse  ratio  implies  a  Gaussian  channel.  Channels  that  have 
nearly  all  the  received  signal  energy  in  the  diffuse  component,  i.e.,  direct-to-diffuse  ratio 
less  than  one,  have  strong  fading.  A  direct-to-diffuse  ratio  of  zero  implies  a  Rayleigh 
channel.  For  direct-to-diffuse  ratios  between  these  extremes  the  channel  experiences 
Ricean  fading.  This  thesis  examines  the  performance  in  each  of  these  broad  categories. 

In  addition  to  these  broad  categories,  the  time  varying  nature  of  the  channel  can  be 
described  as  slow  or  fast.  In  this  thesis  the  channel  properties  are  assumed  to  be  constant 
over  the  duration  of  a  hop  and,  therefore,  slowly  varying.  Further,  the  channel  may 
introduce  some  signal  distortion  arising  from  the  treatment  of  sinusoids  comprising  the 
signal  set  within  a  hop  differently.  This  distortion  is  characteristic  of  frequency  selective 
channels.  However,  it  is  reasonable  to  assume  that  the  signal  sinusoids  experience  the 
same  multi-path  effects.  This  is  the  case  in  frequency  nonselective  channels.  One 
measure  of  this  phenomenon  is  the  coherence  bandwidth  of  the  signal.  The  coherence 
bandwidth  is  the  frequency  range  over  which  the  signal  frequencies  pass  through  without 
distortion.  This  can  be  summarized  mathematically  as 

<40,  =  4-  w 

*  m 
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where  (A/)c  is  the  coherence  bandwidth  and  T„,  is  the  multi-path  spread  of  the  channel. 

A  frequency  nonselective  cannel  displays  a  coherence  bandwidth  that  is  larger  than  the 
signal  bandwidth.  The  rate  of  fading  is  related  to 

M'-i  (5) 

The  coherence  time  is  (A/)c  and  Bd  is  the  Doppler  spread  of  the  channel.  Slowly 
fading  channels  display  a  large  coherence  time  or,  conversely,  a  small  Doppler  spread. 
These  descriptors  are  discussed  further  in  the  system  description. [Ref.  2]  The  use  of 
diversity  to  mitigate  the  multi-path  effects  for  conventional  MFSK  has  been  widely 
investigated  [Refs.  2,  3,  4]. 

When  analyzing  performance  it  is  important  to  distinguish  between  a  constant  energy 
per  hop  system  and  a  constant  energy  per  symbol  system.  As  diversity  increases,  the  total 
symbol  energy  in  a  constant  energy  per  hop  system  increases,  while  the  symbol  rate 
decreases.  However,  in  the  constant  energy  per  symbol  system,  increasing  diversity,  L , 
implies  decreasing  hop  duration,  Th.  Hence,  a  constant  data  rate  is  maintained,  but  the 
energy  per  hop  is  reduced.  Since  many  practical  military  communication  systems  employ 
a  fixed  hop  rate  and  a  variable  data  rate,  the  constant  energy  per  hop  assumption  is  more 
logical  [Ref.  5]. 

C.  CONSTANT  ENERGY  PER  HOP  SYSTEMS 

Previous  work  has  examined  the  performance  of  the  noncoherent  MFSK  receiver  in  a 
Ricean  fading  channel  with  partial-band  noise  interference  and  constant  energy  per 
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symbol  [Ref.  3].  Additionally,  the  performance  of  the  conventional  BFSK  receiver  under 
partial-band  noise  jammin  without  the  multi-path  effect  of  fading  is  analyzed  in  [Ref  6). 
However,  these  investigations  do  not  consider  constant  energy  per  hop  signaling.  The 
performance  of  several  different  types  of  diversity  combining  receivers,  including  the 
conventional  receiver,  utilizing  constant  energy  per  hop  is  simulated  for  Rayleigh  fading 
channels  in  [Ref  7], 

Motivating  this  thesis  is  the  uncertain  degree  of  improved  performance  offered  by 
more  elaborate  receiver  designs,  such  as  the  noise  normalized  receiver,  over  the 
conventional  receiver  for  variable  data  rate  systems.  [Refs.  3,  4]  The  expense  and 
complexity  of  a  more  elaborate  receiver  may  not  be  justified  in  some  circumstances  when 
utilizing  a  constant  energy  per  hop  system. 

The  constant  energy  per  hop  scheme  is  of  practical  value  because  it  allows  the 
potential  for  an  adaptive  signaling  scheme  in  which  the  sender  and  receiver  can  optimize 
the  data  transmission  rate.  In  poor  environments,  diversity  can  be  increased  at  the 
expense  of  lowered  data  rate.  In  favorable  environments,  the  level  of  diversity,  L  ,  can  be 
lowered  to  accommodate  a  higher  data  rate.  These  adaptations  will  ideally  not  require 
any  hardware  modifications  and  will  be  transparent  to  a  channel  observer.  With  this 
adaptive  scheme  in  mind  ,  this  thesis  examines  the  improvement  offered  by  varying  the 
level  of  diversity  possible  when  the  jammer  is  sub-optimal. 
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D.  TONE  INTERFERENCE 

Another  type  of  narrowband  interference  is  tone  interference.  This  can  consist  of  a 
single  interfering  or  multiple  interfering  tones.  In  this  thesis,  a  performance  analysis  for  a 
FFH/BFSK  receiver  without  diversity  over  fading  channel  conditions  similar  to  those 
assumed  for  noise  interference  is  considered.  Since  the  interfering  tone  or  tones  are 
signal-like  in  nature,  they  too  can  suffer  multi-path  effects.  This  analysis  considers  the 
effects  of  fading  on  both  the  signal  and  jammer.  In  a  FFH/BFSK  system  an  intelligent 
jammer  can  potentially  cause  more  degradation  by  splitting  his  available  power  over  both 
the  signal  tones.  The  degree  of  communication  impairment  of  the  single  tone  interference 
per  hop  versus  two  interference  tones  per  hop  strategy  is  also  considered.  Clearly,  the 
greatest  performance  degradation  occurs  when  the  interfering  tones  correspond  exactly  to 
the  various  frequency-hopped  symbol  tones.  Tone  jamming  where  the  tones  do  not 
correspond  exactly  to  the  various  frequency-hopped  symbol  tones  are  not  considered  in 
this  thesis. 
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II.  PARTIAL-BAND  NOISE  INTERFERENCE 


A.  SYSTEM  DESCRIPTION 

The  partial-band  noise  interference,  either  intentional  or  unintentional,  considered  in 
this  thesis  is  modeled  as  additive  Gaussian  noise  and.  when  present,  is  assumed  to  be  in 
all  branches  of  the  MFSK  demodulator  for  any  reception  of  the  dehopped  signal. 

Thermal  noise  and  other  wideband  interferences  which  are  also  assumed  to  corrupt  the 
signal  are  modeled  as  additive  white  Gaussian  noise.  Only  the  signal  is  assumed  to  be 
affected  by  channel  fading.  The  smallest  spacing  between  frequency  hop  slots  is  assumed 
larger  than  the  coherence  bandwidth  of  the  channel,  hence,  each  dehopped  signal  fades 
independently  [Refs.  2.8],  As  discussed  in  the  Introduction,  the  signal  bandwidth  is 
assumed  to  be  much  smaller  than  the  channel  coherence  bandwidth,  and  the  channel 
coherence  time  is  assumed  to  be  much  larger  than  the  hop  duration  or.  equivalently,  the 
hop  rate  is  assumed  to  be  large  compared  to  the  Doppler  spread  of  the  channel.  The  first 
assumption  implies  that  the  channel  is  modeled  as  frequency-nonselective,  while  the 
second  implies  that  the  channel  is  slowly  fading.  The  signal  channel  is  modeled  as  a 
Ricean  fading  channel,  hence,  signal  amplitude  is  a  Ricean  random  variable  [Refs.  2,8]. 
For  Ricean  fading,  the  total  signal  power  consists  of  a  direct  signal  component  and  a 
diffuse  signal  component,  and  the  strength  of  the  fading  channel  is  characterized  by  the 
ratio  of  the  direct  signal  component  power  to  the  diffuse  signal  component  power. 
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The  symbol  rate  is  /?,.  For  MFSK  with  M  order  modulation,  the  corresponding  bit 
rate  is  .Jfj,  =  log  , (M) .  For  L  hops  per  symbol,  the  hop  rate  is  Rh  =  LRS.  The  spread 
spectrum  bandwidth.  W .  is  considered  very  large  compared  to  the  hop  rate. 


B.  PARTIAL-BAND  NOISE  JAMMING  ANALYSIS 
1.  Problem  Development 

The  partial-band  noise  interference  when  present  is  assumed  to  be  in  all  branches 
of  the  MFSK  demodulator  and  affects  each  chip  of  the  dehopped  signal  with  probability  y 
where  y  is  the  fraction  of  the  spread  bandwidth  being  jammed.  Hence,  the  fraction  of  the 

spread  bandwidth  not  being  jammed  is  1  -  y .  If  the  average  power  spectral  density  of 
the  interference  is  Ni/2  over  the  entire  spread  bandwidth,  then  the  power  spectral 
density  of  the  partial-band  interference  when  present  is  vV,  /2y.  The  power  spectral  density 
of  the  thermal  noise  and  other  wideband  interferences  which  are  modeled  as  additive 
white  Gaussian  noise  is  N./2.  Consequently,  total  noise  power  spectral  density  is  Np/2  in 
the  absence  of  partial-band  interference;  otherwise. 


Mg)*® 


(6) 


is  the  total  noise  power  spectral  density  when  narrowband  interference  is  present. 

If  the  equivalent  noise  bandwidth  of  each  detector  branch  in  the  MFSK 
demodulator  is  B  Hz,  then  the  noise  power  received  in  a  given  hop  is  a*  =  g2t  =N05 

with  probability  1-y  when  no  jamming  is  present.  When  jamming  is  present,  the  total 

2  2  2  ( Ni 

noise  power  in  a  given  hop  k  is  ck  =  a,  +  oT  =  +  N0  J  B  with  a  probability  of  y. 
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We  assume  that  each  receiver  hop  slot  has  the  same  noise  equivalent  bandwidth.  The 
noise  equivalent  bandwidth  of  the  receiver  investigated  in  this  thesis  is  B=Rh. 

2.  Probability  of  Bit  Error 

When  partial-band  interference  is  present  the  probability  of  symbol  error  for  a 
MFSK  receiver  is 

L  (  I  \ 

Ps  =  I  .  y'(1  -y)  Ps(i)  (7) 

/=0  V  l 


where  Ps(i )  represents  the  conditional  probability  of  a  symbol  error  given  that  i  of  L  hops 
of  a  symbol  are  jammed.  Since  each  signal  branch  of  the  receiver  is  symmetric  with  the 
other  branches,  we  can  determine  the  probability  of  a  symbol  error  .  Ps .  by  considering 
the  signal  to  be  present  only  in  branch  one  of  the  MFSK  demodulator.  The  outputs  of  the 
other  branches  are  assumed  identical  and  independent  (iid). 

For  orthogonal  MFSK  the  probability  of  a  bit  error  is  related  to  the  symbol  error 
by 


P* 


2{M-  1) 


(8) 


The  energy  per  bit  as  a  function  of  the  symbol  energy  and  the  modulation  order  is 


Es 

log2(A/) 


(9) 
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3.  Probability  of  Symbol  Error  Under  Partial-band  Noise  Interference 


Assuming  the  signal  is  present  in  branch  one  of  the  MFSK  demodulator  allows 
us  to  write  the  probability  of  a  symbol  error  based  on  the  conditional  probability  density 
functions  that  /  of  L  hops  are  jammed. The  conditional  probability  density  function  for  the 

output  of  the  branch  containing  the  signalis  /v,(*i  I/)  where  X,  is  the  random  variable 
that  represents  the  output  of  the  signal  branch.  The  conditional  probability  density 

functions  for  the  non-signal  branches  ar ef\„(xm I i).m  =  2, 3,4..., M  where  the  A'm ’s  are 
the  identically  distributed  random  variables  that  represent  the  output  of  the  branches  that 

do  not  contain  the  signal.  [Ref.  3]  The  conditional  probability  of  symbol  error  is 


PAi)  =  1  -  £/v,(*i  If)  U'/uOn \i)dxm 


tv- i 


dx  i 


(10) 


for  all  m  *  1 . 

Since  the  partial  band  interference  may  or  may  not  be  present  in  a  hop.  we  must 
be  able  to  differentiate  between  the  two  possibilities.  Let  the  subscript  n=  1,  2  denote  that 
hop  k  of  a  symbol  has  interference  and  has  no  interference,  respectively.  The  diversity 
summer  acts  to  add  together  all  L  independent  hops  in  its  branch.  Of  those  hops  i  of  them 
are  jammed;  hence,  we  can  express  the  random  variable  at  the  output  of  the  diversity 
summer  ( 1 0)  as 


/  <  / 

A’m  =  S  A'm*„  =  £  Xmk,  +  X  A  m*-  »»  =  1.2.  3. ...A/ 

*=l  A=l  k—l.—i 


(ID 
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a.  Probability  Density  Function  for  the  Decision  Variable  Xm 

The  probability  density  function  of  the  independent,  identically  distributed 
(iid)  random  variables  X,* ,  m=  2,3,. ..,M  that  represent  the  demodulator  branch  outputs 
not  containing  the  signal  for  hop  k  of  a  symbol  is 


/v^M  =  77exp 
2  at 


'  \ 

Xmk 

<  2&1  J 


u(xmk )  m  =  2,3,...,M 


(12) 


where  w(«)is  the  unit  step  function.  [Ref.  9] 
The  Laplace  transform  of  (12)  is 


1 


2ajs  + 1 


(13) 


Since  each  hop  is  independent,  from  ( 1 1 ),  we  can  express  the  Laplace  transform  of  the 
conditional  density  function,  f\m(xm  I i)  ,  as 

^.(ilO  =  [F.v»,(jlo]'x[^-J(ilo]''"  (14) 


Direct  inversion  of  (14)  leads  to  an  infinite  series  of  confluent  hypergeometric  functions, 
but  this  proved  difficult  to  program  and  slow  to  execute.  To  evaluate  (34)  we  first  invert 
the  individual  portions  consisting  of  either  all  hops  jammed  or  all  hops  not  jammed.  The 
conditional  density  function  for  the  decision  variable  is  then  obtained  by  convolution. 

The  Laplace  transform  pairs  are 


and 


(15) 
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fx^(s\i)]1 


2  olJ  (L-i-  1) 


— — —  exp  —  u(x ) 
—  /  —  1 )!  7<tT, 


where  ct*::  and  ak:2  represent  the  noise  power  in  a  hop  experiencing  interference  and  a  hop 

experiencing  only  thermal  noise,  respectively.  To  recover fxm{xm)  ,  the  two  m-Erlang 
random  variables  above  are  convolved  .  Hence,  the  inverse  Laplace  transform  of  (14)  is 

obtained  from  a  convolution  of  the  individual  right  hand  sides  (RHS)  of  (15)  and  (16)  as 


A.<*.io=JMrJrJ  IrVI  r~ 17  e*p  I  rr  a" -  f-  W exp  1  I  ** ( 1 7) 

V2cA)  /  \2o ft  J  v  *  )•  V 2(J^|  /  \L  t  1 ) .  \  2c  k2 

Applying  the  binomial  theorem  to  the  integrand  of  ( 1 7)  and  defining 


(xm  - 1) 
2cl 


2c\x)  [2 c\2)  {L-i-  !)!(/-  1)! 


we  can  express  ( 1 7)  as 


fxm{xm\i)  =  \  C"  Z 


f  L-i- l 


xkm(-t)L  ,_l"*exp 


_j _ ]_ 

'  2(Jj^  ■ 


dt  (19) 


Interchanging  the  order  of  integration  and  summation  and  integrating  term  by  term,  we 
get  [Ref.  10,  equation  3.351.1] 


/uxJo^'z  ( ■ L-'-' 1  (20) 

*=0  V  *  7  V°-Py  /=0  (k  +  i-  1  - AKa-Br 


where  for  brevity 
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(22) 


P  = 


2cr*2 


The  special  cases  of  all  hops  jammed  (i=L)  and  no  hops  jammed  (;-0)  are 

obtained  directly  from  (15)  and  (16)  respectively,  since  the  analysis  leading  to  (20)  is  not 

valid  in  theses  cases.  The  remaining  task  is  to  evaluate  the  cumulative  distribution 

function  of  fXm(xm\i)  that  is  required  for  computing  the  probability  of  symbol  error,  (10). 

The  computation  required  for  the  non-signal  branch  contribution  to  the  probability  of 
symbol  error  is  fxm(xm).  There  are  two  avenues  to  compute  this,  the  first  based  on 

direct  integration  or,  equivalently,  the  second  based  on  an  inverse  Laplace  transform  of 
F.\m(s\i)/s. 

Directly  integrating  (20),  valid  when  0  <  i  <  L,  we  have 


(a')p 


L-i- 1  ( 

I  H)‘ 


L-i-l 

k 


x(k  +  /-!)!  {L-  ‘  1  k)[  -1-j  -  exp(-p*m/  I  * - *= - — - - 

(a  -  P)*+l  Ip4"-1  0  (L-i-\-n-p)\^ 


(L-i-  1)! 


k+i-) 

j  - 

/= 0  (k  +  i-  1  -l)(a  -  P)/+l  (a 


L-2-m 


-cxp(-axm)  £ 


L-2-l-q 

Xm 


q-o  (L  -  2-  I  -  q)\a^' 


(23) 


This  equation  is  simplified  by  expanding  the  binomial  coefficient  and  canceling  like 
terms  to  obtain 
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p-s  ,  1  ^  (*  +  i-l)!  a'p* 

J0  fxm(xm\i)dxm  =  — —  Zj  - 77 - ( - - -  . 

Jo  O'1)?  *=o  k\  (-!)  (P~a) 


L~l~  1 

l-exp(-jkm)  2 


(P*«)' 


P=o  P'- 


nL-i  k+i- 1 

P  I 


(L-2-  /)!(-!) 


k-l-\ 


(L-i-  \  - k)\  /=0  (*  +  /-  l  -/)!(p-a)m+la/-'-,-/ 


l-exp(W2'^ 

<?=o  P- 


(24) 


when  0  <  i  <  L.  For  the  case  of  all  hops  jammed,  i  =  I,  the  cumulative  density  function  is 


£Vw»*-E-(^r 


y^- 1  y 

exp  (  — —  \dx 


(I-D!  r  2a2 


*i 


=  1  -  exp 


x  Me1  i 


(  \  * 


.y 

.2a*,  J*=o  k\  ^2a*,  j 


(25) 


and  for  the  case  of  no  hops  jammed,  i  =  0,  it  is 


l'h.(x)dx  =  £■ 


2a^  7 


ri-l 


x  exp^M* 


(1-1)! 


2a  *% 


4=J 


=  1  -  exp • 


*  1 


2a*2  J  *=o  \2o*: 


(26) 


A  final  useful  form  is  obtained  by  making  the  linear  transformation  of  the  random 

variable  Xm  .  Recognizing  that  a  =  and  P  =  then 

2a]  2o; 


a.  _  2<Jf  _  £2  _  N0 
a  1 


P  -L  0?  yvr 


(27) 


2a; 


Then  making  the  change  of  variables  y  =  p*m  in  (24)  we  get 
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Fy(yl/)  = 


1 

0-1)! 


(k~  - 1)! 
n\ 


(NtY 

(-1)  \NT~  1)*+' 


L-i-\-k 

1  -exp(-j')  2 


p= o 


(L -2 -/)!(-!  )*~M 

- o?  h  (k+i-\-t)\(NT-\)M 


f  \ 1 

(  y  Y&  UJ 


9! 


(28) 


which  is  useful  for  numerical  integration  purposes.  Applying  the  same  substitution  in  (25) 
and  (26),  we  get 


and 


„a,I  l\  -  1 


£-1 


fY(ylO)  =  1  -  exp 


y_  y  i_(y_)k 
Nt  £0  k\  W 


(29) 


(30) 


for  the  case  of  all  hops  jammed  and  no  hops  jammed,  respectively. 

Previous  efforts  to  determine  the  cumulative  distribution  for  the  non-signal 
branches,  the  second  technique,  is  based  on  a  numerical  inverse  Laplace  transform  of  the 
function  j  x  F,\m(s\i)  [Refs.  3, 4].  Despite  the  finite  sums  of  exponential-like  terms  in 
(28).  this  second  method  proved  more  accurate  for  large  values  of  L  . 

b.  Probability  Density  Function  for  the  Decision  Variable  X, 

Continuing  with  the  assumption  that  the  signal  is  present  in  branch  one  of  the 
MFSK  demodulator  accompanied  by  Gaussian  noise  ,  the  random  variable  X]  is  described 


17 


by  the  conditional  density  function  /v,  (xi  I/) .  Given  a  signal  amplitude  of  Jl  ak  and 


before  diversity  combining,  the  probability  density  function  of  X,  k  is 


_1 
2  o 


2  exp 


(31) 


where  /<)(•)  is  the  modified  Bessel  function  of  the  first  kind  and  order  zero.  The  channel 


is  assumed  to  have  Ricean  fading;  therefore,  ak  is  a  Ricean  random  variable  with  a 


probability  density  function  given  by 


fAt(ak)=- 1  exp 


(32) 


where  a2  is  the  average  power  of  the  direct  signal  component  and  2a2  is  the  average 
power  of  the  diffuse  signal  component  [Refs.  2,  11].  The  total  average  signal  power  of  a 

hop  it  of  a  symbol  is  a2  +  2a2  and  is  assumed  to  remain  constant  from  hop  to  hop. 

The  random  variable  X\k  when  hop  k  of  a  symbol  has  interference  and  no 
interference  is  denoted  by  Xu,  and  X\k2  respectively.  The  conditioning  is  removed  by 

integrating  the  product  of  (31)  and  (32)  with  respect  to  ak  from  zero  to  infinity  with  the 
result 


/v,»„(*u„)  = 


2i°l 


+  2  a2 


exp 


_]_  in,  +  2oc~ 

2  V  a*„  +  2c2  J 


I  o 


^2xH„a2 

o;  +2a: 
*>» 


(33) 
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for  x\k„  £  0  where  a*,  implies  a  jammed  hop  and  cs2k2  implies  no  jamming.  Since  each 
hop  is  independent,  from  (1 1 )  we  see  that  the  Laplace  transform  of  the  conditional 

density  function  for  the  signal  branch,/v,(*i  I/)  ,  is 

=  [F.v.Csl/)]'  •  (34) 


The  Laplace  transform  of  the  density  function  describing  a  single  hop,  /*,*(*  u)  ,  is 

'-2a2p*y 


FXu  (^)  =  exp 

(s  +  PO 


s  +  P  *„ 


(35) 


where 


(36) 


When(35)  to  the  c„  power  and  the  inverse  Laplace  transform  is  taken,  the  result  is 
equivalent  to  taking  the  c„  fold  convolution,  denoted  by  ®c„  .  Letting  c,  =  i  and 
c2  =  L  -  i,  so  that  elements  subscripted  with  a  1  correspond  to  jammed  hops  while 
elements  subscripted  with  a  2  correspond  to  hops  experiencing  only  thermal  noise,  then 


o  (c„-l)/2 

(2c„cr) 


xexp  [-p*„(x u„  +  2c„a2)]  x  1Ch_x  {2^k„j2cna2xik„  )  w(xu„)  (37) 

where  lCn- 1  represents  the  modified  Bessel  function  of  integer  order  c„-l .  [Refs.  3,  1 1  ] 
Unlike  the  case  of  the  non-signal  branch  analysis,  no  analytic  solution  for  the  final 
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convolution  for  the  decision  variable  Xt  has  been  found.  Symbolically,  the  probability 
density  function  for  X,  at  the  output  of  the  diversity  summer  is  described  by 

/v,(xi  1 0  =  [fXui  (xu,  li)]®'  <E>  [/v„2(xu: \i)]m~‘)  (38) 


C.  NUMERICAL  EVALUATION  OF  PARTIAL-BAND  NOISE  JAMMING 

For  levels  of  diversity  up  to  L=  20,  the  conditional  probability  of  symbol  error  is  First 
computed  based  on  the  signal  energy-to-thermal  noise  power  spectral  density  ratio, 
signal-to-jammer  power  spectral  density  ratio,  signal  fading  direct-to-diffuse  ratio,  and 
fraction  of  spread  spectrum  bandwidth  jammed.  These  input  parameters  form  the  basis 
for  the  probability  of  symbol  error  which  is  converted  to  probability  of  bit  error  using 
(8).  To  summarize 

_  a*  direct  signal  power/hop 

S*  =  ~j  = - : - z -  (39) 

a  f  noise  power/hop 


diffuse  signal  power/hop 
noise  power/hop 


(40) 


Eb _ average  energy  per  bit _ 

N0  thermal  noise  power  spectral  density 


(41) 


Eh _ average  energy  per  hop 

Nj  interference  power  spectral  density 


(42) 


and  0.0  <  y  <  1 .0  . 


(43) 
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The  computation  of  the  probability  of  symbol  error  is  achieved  for  each  L  by 
weighting  and  summing  the  probabilities  of  symbol  error  conditioned  on  /  of  L  hops 
jammed.  A  numeric  integration  is  performed  over  the  product  of  the  signal  path 
contribution  and  the  non-signal  path  contribution  (10). 

1.  Non-Signal  Path  Contribution 

The  cumulative  distribution  function  is  computed  using  the  analytic  expression  of 
finite  summations  given  by  (28)  for  0  <  /  <  L  for  small  L.  However,  it  is  computationally 
faster  and  results  in  less  round  off  error  for  large  values  of  L  if  a  numerical  inverse 

Laplace  transform  is  performed  on  ~  x  FXm(s\i).  The  special  cases  of  all  hops  jammed 
(i-L)  or  no  hops  jammed  (r=0)  can  be  computed  directly  from  (25)  or  (26),  respectively. 

The  convergence  of  the  inverse  Laplace  transform  is  accelerated  by  the  Euler 

transformation  [Ref.  12]  usually  taking  about  60  terms  to  reach  relative  errors  on  the 

order  of  1  O'9.  However,  certain  combinations  of  input  parameters  results  in  slower 

convergence.  The  inverse  Laplace  transform  algorithm  is  limited  to  not  more  than  1000 

iterations. 

2.  Signal  Path  Contribution 

For  the  special  cases  either  of  all  hops  jammed  or  no  hops  jammed,  the  required 

probability  function  /*,(*!  I/)  is  given  by  (37).  In  the  case  of  all  hops  jammed  n=l, 
c,  =  Z,,  and  c2  =0,  and  in  the  case  of  no  hops  jammed  n=2.  c,  =  0.  and  c;  -L.  In  both  of 

these  cases,  (10)  can  be  solved  analytically,  but  the  results  are  so  complex  that  numerical 
evaluation  is  easier  and  more  straightforward.  When  i  *  0  and  i  *■  Z, .(38)  is  evaluated  by 
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the  numerical  inversion  of  (34)  after  which  the  probability  of  symbol  error,  (10),  is 
evaluated  numerically. 

The  partial-band  jamming  fraction  that  yields  the  worst-case  performance  for  the 

conventional  receiver  is  obtained  experimentally  by  computing  the  probability  of  bit  error 

Eh  Eh 

as  a  function  of  y  for  fixed  values  of  —  and  — .  To  cover  the  broad  range  of  channel 

N0  Nj 

fading  severity,  results  are  obtained  for  several  values  of  the  direct-to-diffuse  ratio.  For 
weakly  fading  channels,  exhibiting  a  strong  direct  signal  component,  a  direct-to-diffuse 
ratio  of  100  is  used.  This  is  essentially  a  Gaussian  channel.  For  strongly  fading  channels, 
considered  as  nearly  Rayleigh  channels,  a  direct-to-diffuse  ratio  of  one  is  used.  Typical  of 
a  moderate  Ricean  fading  channel,  a  direct-to-diffuse  ratio  of  ten  is  used  . 
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D.  PARTIAL-BAND  NOISE  JAMMING  RESULTS 


To  determine  the  worst  case  partial-band  jamming  performance  of  the  FFH/MFSK 
receiver  the  probability  of  bit  error  is  computed  as  a  function  of  y.  The  process  is  repeated 
for  levels  of  diversity  up  to  twenty,  which  represents  the  upper  end  for  practical  systems 
in  today's  technology.  Figure  2  through  Figure  9  display  the  probability  of  bit  error  for 
various  levels  of  signal-to-noise  ratio,  total  jammer  power,  and  modulation  order.  The 
figures  are  calculated  with  the  assumption  that  the  signal  direct-to-diffuse  ratio  (Rb)  is 
constant  during  a  hop  duration.  Figure  2  demonstrates  the  performance  for  moderate 
signal-to-noise  ratio  with  a  near  Gaussian  channel.  The  signal-to-thermal  noise  ratio, 
Eh/N0 ,  is  the  signal  energy  contained  in  one  hop  for  a  conventional  BFSK  receiver.  The 
table  below  summarizes  the  cases  displayed  in  Figure  2  through  Figure  9. 


TABLE  1.  PARAMETERS  FOR  WORST  CASE  ANALYSIS 


FIGURE 

Eh/NcdB 

Et/NjdB 

R,,  /  Ricean  Fading 

Figure  2 

13.35 

3 

100  /  Very  Weak 

Figure  3 

13.35 

3 

1 0  /  Moderate 

Figure  4 

13.35 

10 

100 /Very  Weak 

Figure  5 

13.35 

10 

1 0  /  Moderate 

Figure  6 

16 

3 

1 0  /  Moderate 

Figure  7 

16 

10 

1 0  /  Moderate 

Figure  8 

20 
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Figure  2.  Worst  Case  Partial-band  Jamming 
Eh/N0=13.35dB,  Eh/A[/=3dB,  R^lOO 

The  surprising  result  is  that  despite  increased  levels  of  diversity,  the  performance 
does  not  greatly  improve.  There  will  be  more  to  say  about  this  later,  but  the  same  trend  is 
observed  in  Figure  3  where  there  is  Ricean  fading  of  the  signal . 


2' 


0.300 

0.200 

0.100 

0.050 

0.030 

1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  19  20 

Number  of  Hops 

Figure  3.  Worst  Case  Partial-band  Jamming 
Eh/N=13.35dB,  E„/A>3dB,  R^lO 

Another,  surprising  result  is  that  increasing  the  modulation  order  to  four  or  eight 
does  not  provide  much  improvement.  Figure  2  and  Figure  3  display  pessimistic 
performance  from  the  communicator's  viewpoint.  However,  they  assume  a  high  level  of 
jammer  power.  We  next  investigate  the  performance  when  the  jammer  power  is  just  one 
tenth  of  the  communicator's  power  when  L=\  (  Figure  4  and  Figure  5) .  Since  we  are 
assuming  1000  hop  slots,  this  is  equivalent  to  assuming  a  signal-to-jammer  power  ratio  of 
-20dB  for  BFSK  with  no  frequency-hopping 
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Figure  4.  Worst  Case  Partial-band  Jamming 
Eh/No=13.35dB,Eh/A(/=10dB,Rb=100 


Number  of  Hops 


Figure  5.  Worst  Case  Partial-band  Jamming 
E|/N0=13.35dB,  Eh/A(/=10dB,  1^=10 


A  comparison  of  Figure  4  with  Figure  5  illustrate  the  influences  of  Ricean  fading 


on  the  probability  of  bit  error  when  the  communicator  enjoys  a  significant  power 
advantage  over  the  jammer.  Next  we  consider  a  situation  where  the  communicator  has  a 
larger  initial  signal-to-thermal  noise  ratio,  i.e.  16dB,  to  see  if  at  higher  levels  the  jammer 
could  mitigate  the  effects  of  diversity  with  partial-band  noise  jamming.  Figure  6  and 
Figure  7  display  these  results. 
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Figure  6.  Worst  Case  Partial-band  Jamming 
Eh/N0=  1 6dB,  Eh/iV/=3dB,  1 0 


27 


Number  of  Hops 

Figure  7.  Worst  Case  Partial-band  Noise  Jamming 
Eh/N^lbdB,  Eh/A//=10dB,  R^IO 

To  examine  the  extremes  of  channel  fading  on  the  received  signal,  the  case  of 
near  Rayleigh  fading  is  examined  with  a  direct-to-diffuse  ratio  of  one.  Since  fading  is 
severe,  we  would  expect  poor  performance  even  without  jamming,  but  with  some 
improvement  added  by  the  increased  diversity.  The  signal  energy  per  hop  is  examined  at 
20  dB  (Figure  8)  and  at  lOdB  (Figure  9). 
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Figure  8.  Worst  Case  Partial-band  Jamming 
E„/No=20dB,  Eh/A0=3dB,  R^l 
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Figure  9.  Worst  Case  Partial-band  Jamming 
Eh/No=10dB,  Eh/iV/=0dB,  Rt=l 


In  both  Figure  8  and  Figure  9  the  improvement  offered  by  diversity  initially 
lowers  the  probability  of  bit  error,  but  not  to  the  extent  that  it  can  overcome  the  errors 
introduced  by  partial-band  jamming.  It  should  be  noted  also  that  the  high 
signal-to-thermal  noise  ratio  is  received  after  passing  through  the  near  Rayleigh  channel. 
Therefore,  it  has  a  diffuse  signal  component  equal  to  the  direct  signal  component.  This 
represents  a  very  poor  channel  for  the  communicator  that  is  further  degraded  by 
partial-band  noise  jamming. 

The  conclusion  to  be  drawn  from  the  preceding  series  of  performance  plots  is  that 
for  a  wide  range  of  signal-to-thermal  noise  level,  jammer  power,  and  channel  conditions, 
the  increase  in  diversity  in  the  FFH/MFSK  receiver  alone  is  not  enough  to  overcome  the 
effects  of  pariial-band  noise  jamming  acting  in  concert  with  multipath  effects.  The 
critical  parameter  then  is  the  fraction  of  the  spread  spectrum  bandwidth,  y,  over  which  the 
jammer  power  is  spread.  The  Figures  10  through  17  demonstrate  the  fraction  of  jammed 
bandwidth  that  provides  the  worst  case  probability  of  bit  error. 
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Figure  10.  Fraction  of  Jammed  Bandwidth 
Eh/N0=13.35dB,  E/A>=3dB,  R„=100 
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Figure  1 1 .  Fraction  of  Bandwidth  Jammed 
Eh/N0=13.35dB,  Eh/A0=3dB,  1^=10 
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Figure  12.  Fraction  of  Bandwidth  Jammed 
Eh/N0=13.35dB,  Eh/N/=10dB,  1^=100 
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Figure  13.  Fraction  of  Bandwidth  Jammed 
E„/N0=13.35dB,  Eh/Nj=\0dB,  1^=10 
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The  worst  case  values  y  are  determined  by  computing  the  probability  of  bit  error 
as  a  function  of  y.  Figure  10  through  Figure  17  show  that ,  in  general  y,  depends  on  the 
hop  energy-to-thermal  noise  ratio  (Eh/N0 ),  the  hop  energy-to-jammer  power  spectral 
density  (E  JNj),  the  level  of  diversity  (L),  the  modulation  order  (M),  and  the  channel 
fading  direct-to-diffuse  ratio.  Further,  the  shape  of  the  worst  case  y  is  generally  the 
same,  showing  an  initial  drop  as  L  increases  with  a  gradually  lowering  slope.  It  is 
important  to  note  that  the  lowest  levels  of  y  are  still  above  10'3  which  represent  only  one 
jammed  hop  slot  when  there  are  1000  hops  in  the  frequency-hopping  system. 

That  the  increase  in  diversity  is  unable  to  mitigate  the  partial  band  jamming  is 
also  indicated  by  looking  at  the  trend  of  performance  for  fixed  values  of  y.  Figure  10 
through  Figure  17  show  that  performance  improves  linearly  as  diversity  increases.  This 
is  not  unexpected  since  the  constant  energy  per  hop  system  increases  the  total  energy  per 
bit  as  diversity  increases,  however,  the  change  in  slope  of  the  performance  curves  is 
indicative  of  the  effect  contributing  to  the  worst  case  analysis.  Figure  1 8  through 
Figure  23  display  the  performance  for  fixed  values  of  y. 
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Figure  18.  Performance  for  fixed  values  of  gamma 
Eh/N0=13.35dB,  Eh/A^10dB,  R„=10 


Figure  18  displays  the  performance  in  a  moderate  Ricean  fading  channel  where 
communicator  enjoys  a  fair  signal-to-jammer  power  advantage.  The  performance 
improves  linearly  as  diversity  increases.  However,  the  slope  of  the  performance  curve  is 
significantly  reduced  for  decreasing  values  of  y.  Further,  we  can  see  that  for  a  given  L 
there  corresponds  a  worst  case  value  for  y  which  becomes  less  sensitive  as  diversity 
increases.  Figure  19  displays  the  performance  when  the  signal  enjoys  less  power 
advantage  relative  to  the  jammer.  The  trend  is  that  for  lower  signal-to-jammer  power 
ratio  performance  improves  more  slowly. 
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Figure  19.  Performance  for  fixed  values  of  gamma 
Eh/N0=l 3.35dB,  Eh/A(/=3dB,  R„=10 


The  worst  case  fraction  of  bandwidth  jammed  is  influenced  by  channel  fading 
(Figure  10  through  Figure  17).  Likewise,  the  slope  of  the  performance  curves  is 
influenced  by  channel  fading  .  As  the  fading  becomes  more  severe,  i.e.,  near  Rayleigh, 
the  performance  improves  more  slowly  as  demonstrated  in  Figure  21  through  Figure  23. 


37 


III.  MULTI-TONE  INTERFERENCE 


A.  SYSTEM  DESCRIPTION 

In  addition  to  partial-band  noise  jamming,  spread  spectrum  communications  can 
experience  tonal  or  narrowband  interference  in  one  or  more  hop  slots.  We  can  similarly 
expect  that  as  spread  spectrum  bandwidth  grows  the  likelihood  of  tonal  interference  will 
also  grow.  This  interference  may  arise  from  hostile  sources  as  well  as  non-hostile  ground 
communications,  satellite  transponders,  and  radars.  The  multi-tone  interference  problem 
is  similar  to  the  partial-band  noise  problem. 

This  chapter  analyzes  the  probability  of  bit  error  arising  from  an  "intelligent  jammer" 
who  knows  the  spread  spectrum  bandwidth,  W,  and  the  modulation  order.  This 
intelligent  jammer  is  able  to  place  a  single  interfering  tone  in  one  or  more  hop  slots,  but 
has  no  knowledge  of  the  pseudonoise  (PN)  sequence  driving  the  hopping  pattern.  This 
represents  a  very  intelligent  jammer.  We  then  relax  this  situation  and  compare  it  with  the 
performance  of  the  jammer  who  places  an  interfering  tone  in  none,  one,  or  both  of  the 
symbol  tone  locations  in  a  particular  hop  slot. 

It  is  likely  that  the  interfering  tones  may  experience  multi-path  effects  in  a  manner 
similar  to  the  signal  tone.  However,  since  the  interfering  tones  and  signal  tone  can  arrive 
from  different  paths,  the  interfering  tones  will  in  general  arrive  at  the  receiver  with  a 
different  direct-to-diffuse  power  ratio  from  that  of  the  signal.  In  this  thesis,  the 
performance  when  both  the  signal  and  the  jammer  experience  fading  is  evaluated. 
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B.  MULTI-TONE  INTERFERENCE  ANALYSIS 


1.  Problem  Development 

First  we  consider  the  case  where  the  jammer  places  at  most  one  tone  in  each  fast 
frequency-hop  slot.  Only  the  worst  case  situation  where  the  jamming  tones  exactly 
coincide  with  frequency-hopped  symbol  tones  is  considered.  The  situation  may  be 
visualized  with  the  aid  Figure  24. 

FFH/BFSK 
N  slots 


Figure  24.  Multi-tone  interference  in  FFH/BFSK 
The  FFH/BFSK  receiver  is  a  special  case  of  the  FFH/MFSK  receiver  described 

earlier.  The  orthogonal  design  is  selected  such  that  ■'/} b  is  the  system  bit  rate  and  each 

signal  tone  is  placed  b  apart.  In  Figure  24  the  bandwidth  of  a  single  hop.  B,  is  2-'A\ . 

The  spread  spectrum  encompasses  W  Hz  and  is  divided  into  N  hop  slots.  Although 

numbered  consecutively  for  illustration,  the  pseudonoise  modulation  will  create  a 
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random-like  sequence  to  the  hopping  pattern.  In  this  thesis,  N=1000  is  choosen  as  a 
representative  number  of  hop  slots.  Further,  the  number  of  interfering  tones,  q,  must  be 
an  integer 

1  <  q  <  N  (44) 


Since  the  jammer  spreads  his  available  power,  J,  equally  among  q  tones,  each 
interfering  tone  has  a  power  of  J^  =  J/q.  Earlier  work  focuses  on  the  simplifying  case 
of  no  thermal  noise  [Ref.  13].  In  the  no  thermal  noise  case  an  error  can  only  occur  if  the 
hop  is  jammed  and  the  the  interfering  tone  power  is  larger  than  the  signal  power.  In  this 
instance  the  jammer  can  make  most  efficient  use  of  his  power  by  choosing  the  number  of 
tones  to  jam  as 


q  =  INTG 


(45) 


where  Pc  is  the  carrier  power  of  the  desired  symbol  and  INTG[  ]  represents  the  greatest 
integer  less  than  the  argument.  [Ref.  13,  pp. 596-598]  In  this  thesis,  the  probability  of  bit 
error  is  obtained  using  fixed  values  of  q  to  determine  a  worst  case  scenario  which 
includes  the  effects  of  both  signal  and  jammer  fading  as  well  as  thermal  noise.  The 
approach  to  this  problem  is  to  first  detemine  the  probability  of  bit  error  for  a  conventional 
BFSK  system  without  frequency-hopping  and  then  generalize  those  results  to  the 
FFH/BFSK  system. 
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2.  BFSK  Analysis  with  Single  Tone  Interference 


For  the  binary  case  we  begin  with  the  familiar  expression  for  the  probability  of 
bit  error  for  the  noncoherent  FSK  receiver 

Pb  =Ps ,  *Pr[P2  +  x /»,.[/?  i  >/?2l$2]  (46) 

where  R,  and  R2  represent  the  demodulator  outputs  of  the  branch  detecting  signal  5,  and 
s2 ,  respectively.  Making  the  assumption  of  equally  likely  signaling,  we  have 

Ps^PS2  =  \  (47) 

and 

/>6  =  ±[fM*2>/G) +  />,(*!  >*2)]  (48) 

The  received  energy  consists  of  signal  energy,  thermal  noise,  and,  if  jammed,  the 
interference  signal 

r(t)  =  s,(t)  +  n(t)  +  sj(t),  (49) 

where  s ,  is  the  ith  signal  tone,  s,  is  the  interference  tone,  and  n(t)  is  additive  white 
Gaussian  noise  with  power  spectral  density  N0  /  2.  Since  the  receiver  structure  is 
symmetric,  we  can  simplify  the  analysis  problem  by  considering  the  case  where  branch 
two  is  jammed. 

First  we  examine  the  case  where  the  signal  and  interference  occupy  different 
frequencies.  Then 

Pr(error\si)  =  Pr(Ri  >  R\ Isl)  (50) 
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We  express  the  received  waveform  as 


r(t)  =  J2  ac cos  (co  1  /  +  0)  +  v^2  ajcos  (©2  t  +  Bj)  +  n(t) 


(51) 


where  ac  and  aj  are  modeled  as  Ricean  random  variables.  As  previously  discussed  the 
received  amplitudes  and  phases  of  each  tone  is  a  random  variable  due  to  multipath 
effects.  For  the  noncoherent  detector  design,  the  output  of  each  branch  can  be  expressed 
as  a  non-central  Chi  square  random  variable  [Ref.  2].  For  the  signal  branch 


fRi(r\\ac,si)  = 


2  a2 


exp  <  — 


(n+2  a?Tl) 


2a2 


/o 


acTb 


L  a 


J2n 


u(n) 


(52) 


and  for  branch  two  containing  the  interference  tone 


fR1(r2\aJ,sl)  = 


2a 


2  exP  1 


(' -2+2ajli ) 


2ct2 


/o 


ajTb 


J2r2 


u(r2)  (53) 


where  a2  =  N0Th.  The  probability  of  error  (50)  is 

Pr(R2>  R\\ac,ajS\)=  ^R  jR,R1(r\r2\aJ,ac,S\)dr2dr\ 


(54) 


This  represents  the  probability  conditioned  on  the  amplitudes  of  ac  and  a, .  The 
orthogonal  receiver  structure  and  independent  fading  of  the  interference  tone  and  the 
signal  tone  allow  the  joint  probability  density  function  to  be  separated.  Hence.  (54) 
simplifies  to 


Pr(R2  >R\\s])=  j^//?,(ri 


fR2(r2\aj,S\)dr2 


(55) 
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The  evaluation  of  (55)  is  simplified  by  applying  the  quadratic  transformation  to  convert 
(52)  and  (53)  to  Ricean  random  variables  [Ref.  11].  If  r\  =  v?  ,  then  — -  -  2vi .  Using 


this  substitution  in  (52),  we  obtain 


/r,(v,  =  ^TexP - ^2 - |/o 


Similarly  for  the  branch  output  containing  the  interference  tone 


(v?  +  2a?7?)}  |  j2acTbv 


u(y  l)  (56) 


fi,(v2\aj,si)=  ^exp- 
2  CT 


(vi+2  «57») 


- - -  m(v2) 

/TX 


Now  (55)  can  be  expressed  as 


^(^2  >  Filj|)=  J^/r,(V|la0J|)  J^/K2(V2k/,J|)<fr2  ^Vi 

The  inner  integral  equation  has  no  simple  solution,  but  is  represented  symbolically  as 
Marcum's  Q  function  which  is  defined  as  [Ref.  14] 


£?(a,P)=  J?vexp  “(v2  +  a2)  /0[av]rfv 


Hence,  (58)  is  expressed  as 


Pr(v2  >  V|lsi)=  Q^JlajTblo^)  X  li-expl-— -2^— ^  l/0  —  dv]  (60) 


which  can  be  evaluated  to  obtain  [Ref.  1 5] 


Pr(y !  >  v,ls,)  =  i[i  -g(2^.2^)  +e(2^i.2^i) 
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Equation  (58)  represents  the  contribution  to  the  probability  of  bit  error  when  the  signal 
and  jammer  occupy  different  tones  within  a  hop.  It  remains  to  determine  the  contribution 
due  to  a  collocation  of  interfering  and  signal  tones.  That  is 

Pr(e rror  I.S2)  =  Pr(R  1  >  Rj  ^2)  (62) 

where 

/*('.)  =  ^«p{^  ■«(,■,)  (63) 

Since  the  signal  and  the  interference  are  in  the  same  tone  location,  the  output  of 
the  in-phase  and  quadrature  portions  of  branch  two  due  to  signal  and  jamming  tones  is 

J2  Tb  [<arcCOS  (0)  +  flyCOS  (0y)  (64) 

and 

J2  7),[acsin  (0)  +  ajsin  (0y)]  (65) 

respectively.  The  phase  of  the  jammer  is  generally  different  from  that  of  the  signal. 
Squaring  and  summing  to  form  the  branch  output  produces 

^2(0  -  0y)  =  27^[ac  +  +  2  acaj(cos  0  cos  Qj  +  sin  0  sin  0./]  (66) 

which  is  simplified  to 

i?2(8)  =  27’^a^  +a3  +  2ac07Cos5]  (67) 
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Pr(R\  >J?2^2, 6)  -•  J^//?2(nlflc,^',5,52)  x[j^/*,(n  \s2)dr\  dr2  (69) 


Thus,  the  probability  of  error  arising  from  collocation  of  the  jammer  and  the  signal  is 
computed  from  a  probability  density  function  conditional  on  5.  In  this  case  the 
non-signal  branch  is  unaffected  by  the  jamming  tone.  Substituting  (63)  and  (68)  into 


(69),  we  get 


Pr(R\  >^?2 1^2,5)  =  ^xexpi — ~\(al  +a]  +  2aca./cos8) 

2  I  2a 


This  can  be  simplified  into  familiar  terms  by  replacing 


‘2<y2  °'No  N. 


and  similarly  for  the  jammer  power  term 


2  T'b  2  Tb 

a , — -  =>  a /— 
J2  a2  J  N0 


Then  the  conditional  expression  in  (70)  then  becomes 

Pr(R\  >R2\s\,b)  =  exp {y(r*+  Y./  +  2 Jybyj  cos6 ) 


The  nuisance  parameter  5  may  be  removed  by  multiplying  by  its  probability  density 
function  and  integrating  the  conditional  density.  That  is 


Pr(R\  >R2\s\)  =  ~  exp  j  —  (y h  +  yj  +  2  Jyhy j  cos 5 
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where  8  is  modeled  as  a  uniform  random  variable.  Since  the  integral  is  only  over  8,  then 
[Ref.  10] 


Pr(R\  >/?2^i)  =  exp  Y^b+yj)  x  ^  U*  exp  { Jyby j  cos 8)^8  (75) 


which  can  be  evaluated  to  obtain 


Pr(Ri  >R2\si)=^j  xexp  j-i(y6+y^)j  xI0\_Jybyj  ]  (76) 


Therefore,  the  phase  relationship  between  the  signal  and  the  jammer  plays  a  significant 
role  as  evidenced  by  the  modified  Bessel  function  term.  This  makes  physical  sense. 
When  the  phase  is  the  same  or  nearly  so  there  is  reinforcement  of  the  signal  tone  by  the 
interfering  tone,  but  when  the  phase  difference  approaches  1 80  degrees  there  is 
destructive  interference. 

This  allows  us  to  compute  the  probability  of  bit  error  in  the  absence  of  fading  by 
substituting  (61)  and  (76)  into  (48). 


Pb{ac,  aj)  =  i 


+exp 


(77) 


An  equivalent  expression  using  (71 )  and  (72)  is 


Ph(yb,yj)  =  ^  l -Q{jyb,jyj)  +e*p  y(Y6  +  Y./)  h[jy^yj]  (78) 


Using  a  Q  function  property  [Ref.  15,  pp.  396],  we  can  reduce  this  to  an  expression 
containing  just  one  Q  function  as 
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(79) 


Pb(flc,aj)  = 


J 


acajT2b 


Now,  the  amplitudes  of  the  received  interference  tone  and  signal  tone  are  Ricean  random 
variables  where 


fA  ,{aj)  =  exp 
Vj 


r  . 

(  2  ,  ^2^ 

) 

“ 

J 

aj  +  aj 

x/0 

etjctj 

l 

V  2  dj  J 

J 

u(aj) 


(80) 


a. 


fAc(ac)  =  -7  exp  <  - 


dr 


/  9  T  \ 

aj+aj  ' 

ix/0 

acac 

V  2 a]  J 

.  °2e  . 

u(ac) 


Where  in  each  case  the  second  moment  can  be  expressed  as 

a]  =  a]  +  2o]  for  i=  J,  c. 


(81) 


(82) 


We  first  remove  the  conditioning  on  ac  in  (79).  The  technique  is  to  take  each  piece  of 
(79),  multiply  by  (8 1 ),  and  integrate  over  all  possible  values  of  ac.  Hence,  for  the  first 
piece 


j :q[J¥.J¥ 


Ua,  ( ac)dac 


|l 

(  1  ■>  \ 

a c  +  ac 

l 

\  2d]  ) 

Ucdc 

CTc 


(83) 
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olfl 


o2c7i 


+  gx 


l-Q 


olITi 


2  rp2 

ajTb 


olTl 


+  cr^ 


C2cTl  +CT2 


G2cli  +  C2 


Q 


/ 

alfi 

1  a2T2 

u 

a27^  +  a2  ’l 

K^  +  a2  J 

(84) 


The  integral  of  the  product  of  the  second  piece  and  (81)  is 

acacT2b 


2  at 


ac 

x  —  exp  *j  - 

CJ  r 


< a]  +  a 2' 


K 


2a2 


/o 


CtcCLc 


dac  (85) 


Making  the  substitutions 


a  = 


2a2  2a2 


,  B  =  and  T  =j- §,  then 

cr  at 


B2  -  F2  _  +  «cg4  . ,  2a4a^ 

4a  4  a4ac  T^a2  +  a2 


(86) 


at 


jacajtl 


2d  rfc2+<*2 

in  (85),  we  can  use  [Ref.  10,  equation  6.633.4]  to  obtain 


•exp 


-T* 

•  x  /o 

T\aj<xc 

27^a2  +  2  a2 

T'l  ■»  2 

_  T;a;  +  a2  J 

7>2  +  a2 

for  the  second  piece  of  (79).  The  remaining  piece  to  be  determined  is 


fexp{-^(a3+a|)llexp[ 


a2  -f  a2^ 
^  2a2  J 


/o 


acac 


da 


(87) 


(88) 


(89) 
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Using  [Ref.  10,  equation  6.631.4],  we  evaluate  (89)  to  obtain 


7lcr2  +  a; 


exPi 


o)T\  ,  O-lTj 
2°2  2(r^+cJ) 


(90) 


Finally,  combining  (84),  (88),  and  (90)  with  (79),  we  get  the  probability  of  bit  error  when 
the  signal  experiences  Ricean  fading  as 


Pb(aj)  = 


otn 


l{olT2b  +  a2) 


f 

l -Q 

V 

a]fi 

I  a)li 

\a2c 7l  +  o2  *1 

i<y2c7i+<y2  J 

4(7?a?+  a2) 


exp 


a  ffi 

2(7lc?  +  a2) 


exp 


2ct2 


-exp 


n2T2 
aJJ  ft 


}  *  Iv 


ajQicTl 

7fc2+CT2 


(91) 


which  is  now  conditional  on  aj.  Expressing  (86)  in  terms  of 


r  «2^  a2cTh 
^  o2  N„ 


(92) 


c  2a2  T\  2  a2Th 


(93) 


and  the  direct-to-diffuse  power  ratio 
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(94) 


we  get 


Ph(aj)  = 


\b>2 

ft* +  2) 


f 

1-0 

V 


l2aj7i/o2 
\  ^*  +  2 


+ 


1 


2ft* +2) 


0 


2a^/a2  I  2& 


£,*  +  2  ’Ub  +  2 


+ 


1 


2(^6  +  2) 


exp)  - 


C* 


ft* +  2). 


-a]f\\  f  a.jT2h/a2  )  , 

expii^"rexprir^rr/o 


Ta/ct 

^6+2 


(95) 


Using  [Ref.  15] 


2a^/a2  (20r 


£,6  +  2  ’  \  ^6  +  2 


=  1-0 


I  2fr  f2a^/a2 

^6  +  2  ’V  ^6  +  2 


l-> 

(  2a2jT2h/c2+2t;h) 

|  X  /0 

2Kb 

-KrajTh/o 

2 

y  £,6+2  y 

J 

^6  +  2 

L  J 

(96) 


in  (95),  we  obtain 


/>6<fl./)  =  ^ 


1-0 


2^a  2^7^/ct2 


4a  +  2'V  ^6  +  2 


I 


exp- 


-^6 


2(^6  +  2)  '  UU  +  2 


exp  ( -a] T\l2<32 }  +  exp  ■  q2(^^~  '  x  7° 


2  ajThJC^b 
a(^  +  2) 


(97) 
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We  now  employ  a  similar  approach  to  remove  the  conditioning  on  a,.  Integrating 
the  product  of  (80)  and  (97)  over  all  possible  receiver  jammer  amplitudes,  a,,  we  get 


pb  „  1  i-o  - i -  + 

2  l?U(54+5J  +  2)’V(5i  +  5j  +  2  )J 


(4j  +  2)  ..  f-Kt  +  Q)1  .[  2^50 

2(4*  +  4./ +  2)  Ui+4j  +  2j  0  4* +  4-* +  2 


where 


a  )r\  CLjTb 


2  o]Tl  2  a)Tb 


(100) 


R  “I 


(101) 


An  alternate  expression  for  the  probability  of  bit  error  in  the  case  of  both  signal  and 
jammer  fading  in  terms  of  weighted  Q  functions  is 


pb  =  l  i_0  / - S' -  + 

2  nv^+^+2  is**2 ') 


if  J  2 U~  I  2^  )  .[  1  f  2/g^J  ]] 

2  0lu*  +  ^  +  2  V^  +  ^  +  2'J  eXp|  ^  +  ^  +  2j  x/°  ^6  +  ^  +  2 


However,  since  there  is  no  closed  form  solution  for  the  £?  function,  in  practice  it  is  easier 


to  compute  the  performance  using  (98). 
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a.  Special  Cases  of  BFSK  with  Single  Tone  Interference 

The  final  aspect  of  the  analysis  is  to  consider  how  the  derived  result 
simplifies  for  special  cases,  both  as  a  check  on  earlier  work  and  to  see  limiting  cases.  As 
check  on  (98),  when  there  is  no  signal  the  probability  of  bit  error  reduces  to  one  half  as 
expected. 

Next  consider  the  case  of  no  jamming.  In  this  case,  C,}  =  0,  and  the 

probability  of  bit  error  reduces  to 


Pb  = 


x  exp 


(103) 


which  is  the  usual  BFSK  result. 

Another  limiting  case  is  the  performance  when  both  the  signal  and  the 
jammer  experience  Rayleigh  fading.  In  this  case  both  <^b  and  Q  are  zero  and  all  the 
received  energy  is  in  the  respective  diffuse  component.  In  this  case  (98)  reduces  to 

^  +  2 
5*  +  i/  +  2 


(104) 


This  displays  an  inverse  linear  relationship  similar  to  Rayleigh  fading  without  jamming 
[Ref.  2].  Performance  is  improved  if  we  change  the  case  slightly  to  that  where  the 
jammer  suffers  Rayleigh  fading  while  the  signal  experiences  Ricean  fading.  Then 


Pt>  =  ± 


if  i/+2 


2k^  +  4,  +  2 


x  exp 


5*  +  i/  +  2. 


(105) 


Here  any  direct  signal  component  that  reaches  the  receiver  will  serve  to  drive  the 
probability  of  bit  error  lower  in  an  exponential  fashion. 
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Perhaps  the  most  optimistic  performance  from  the  communicator's  point  01 
view  is  when  the  signal  has  no  fading,  but  the  jammer  suffers  Rayleigh  fading.  In  this 
case,  =  £,b  =  0,  and  (98)  reduces  to 


Pb  = 


(106) 


On  the  other  hand  ,  the  most  pessimistic  performance  is  likely  to  occur  when 
the  signal  suffers  Rayleigh  fading,  but  the  jammer  has  no  fading.  In  this  case, 

=  0,  and  (98)  reduces  to 


1- 


1 


&  +2) 


exp 


(107) 


This  performance  will  be  very  poor  unless  there  is  a  relatively  large  amount  of  diffuse 
signal  energy  received. 

3.  FFH/BFSK  Analysis  with  Multi-Tone  Interference 

We  now  turn  our  attention  to  the  application  of  (98)  to  spread  spectrum  and  the 
partial  band  multi-tone  jamming  scenario.  Suppose  N  is  the  number  of  hop  slots  in  the 
FFH/BFSK  system  and  the  jammer  transmits  q  interfering  tones.  When  we  credit  the 
jammer  with  only  placing  at  most  one  tone  per  hop  slot  then 

r=(£)  (108) 

is  the  fraction  of  spread  spectrum  bandwidth  jammed  where  y  is  restricted  to 

i<y<1.0  (109) 

Now, 
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Pi b  =  [error  I  hop  jammed]  +  ^  j  Pr  [error  I  hop  not  jammed] 

where 

/V  [error  I  hop  not  jammed]  =  /^[errorl5FSAT]  =  *  1  exp  j  — —  1  (HO) 

\c,b  +  2J  (.£/>  +  2J 

To  apply  our  earlier  result  for  the  BFSK  receiver  we  must  define  our  jammer  parameters 
on  a  per  hop  basis.  If  y y  is  the  average  total  jammer  power,  then  the  average  total 
jammer  power  per  hop  is 


Y  Y  +  Y  ' 


(111) 


which  implies  and  C,,  are  replaced  by  N^,/q  and  NC,,/q  ,  respectively,  in  (98).  Now  (110) 


is  given  as 


p^{tn)  '-2 


2C, 

^  +  &j  +  2’U»  +  7^  +  2 


^  (;4^2)  [-(;;»  +  fo)  |  ' 

2A^  2(iu  +  £  5, + 2) exp  U*  ■ +  & + 2  J  1 5*  +  &/ + 2 


^  +  2 


(112) 


In  the  case  of  FFH/BFSK  with  one  hop  per  bit,  no  fading,  and  no  thermal  noise 


(N0=4)),  an  analytic  expression  for  the  worst  case  number  of  jamming  tones,  q0 ,  is 
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l<INTG(ffi)=?„<W 


013) 


One  aim  of  this  thesis  is  to  determine  if  this  limited  analytic  result  for  q0  remains  true  for 
FFH/BFSK  over  Ricean  fading  channels  in  the  presence  of  thermal  noise.  This  is 
accomplished  by  comparing  the  performance  based  on  q0  and  some  offset  values  of  q. 

4.  FFH/BFSK  Analysis  Allowing  Two  Interference  Tones 

We  now  extend  the  previous  results  to  a  situation  which  includes  the  possibility 
that  both  receiver  branches  contain  an  interfering  tone.  This  is  a  relaxation  of  the 
intelligent  jammer  scenario  employed  thus  far.  In  this  situation,  the  probability  of  bit 
error  may  be  expressed  as 

Pb  =  /\[no  hops  jammed]  x  Pr  [error  I  no  hops  jammed] 

+/\[one  tone  jammed]  x  Pr  [error  I  one  tone  jammed] 

+/\[two  tones  jammed]  x  Pr  [error  I  two  tones  jammed]  (114) 

The  first  two  conditional  probabilities  were  previously  determined  in  Section  2. 

However,  the  probability  of  those  events  occurring  is  a  sample  without  replacement 
situation.  That  is 

Pr[one  tone  jammed]  =  P,[one  tone  jammed  n  one  tone  not  jammed] 

Pr[0  tone  jammed  ll  tone  not  jammed]  x  Pr[l  tone  jammed  10  tone  not  jammed] 
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(115) 


(  jlVi  _  SzJJ  =  f_L)  (w-i) 

V2 nA  2N-  1 J  \2NA2N-i) 


Similarly 

Pr[ no  tones  jammed]  =  (l  -  (l  -  -) 

_(2N-q\(2N-  1  -q\ 
l  2N  A  2N- 1  J 

And  lastly,  the  probability  that  both  tones  are  jammed  is 

Pr[  both  jammed]  =  Pr[  1  jammed  r\  0  jammed] 

_( q  y  q-l  } 

^2nA2N-  \) 


(116) 


(117) 


All  that  remains  is  to  determine  the  probability  of  error  given  that  both  branches  are 
jammed. 

a.  Extension  of  Single  Interfering  Tone  Results 


Now  that  both  branches  of  the  BFSK  receiver  contain  jamming  tones, 
symmetry  is  restored  .  [Ref.  11]  Hence,  we  can  determine  the  probability  of  bit  error 
supposing  that  signal  s,  is  sent.  An  error  occurs  when  the  output  of  branch  two  is  greater 
than  the  output  of  branch  one  and  (54)  yields 


^[errorltwo  tones  jammed]  =  : 
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Alternatively,  this  may  be  expressed  as 


_  r  i  ,  ( a  iTh  TbJa]+a]  +  2acajcos  5 

Pr [error ! two  tones  jammed]  =  Q  — - - — - 


( -lexn  < 

T\  (a]  +  2a]  +  2acajcos  8  ] 

,  V  1  r\ 

T\aj  Ja]+a]  +  2acajcos  8 

^2;exp 

2  o2 

x  /o 

CT2 

(119) 


The  conditioning  on  the  received  amplitudes  of  the  signal  and  the  jammer  is  removed  by 
multiplying  ( 1 1 9)  with  (80)  and  (81)  and  integrating  over  all  values  of  a,  and  ac .  In 
addition,  the  conditioning  on  the  signal-to-jammer  phase  difference,  8,  must  also  be 
removed. 

In  order  to  make  (1 19)  valid  for  FFH,  we  substitute  (111)  into  (119).  Now, 
symbolically  the  probability  of  error  in  the  case  of  both  tones  jammed,  Ph2 ,  is 


lo 


2\ 

9  ^  _2  ,  2A^  2  ->  I2N  , 

ac  +  1TaJ  +  2  acaj\  —  cos  8 


dacda jdb 


(120) 
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Substituting  our  previous  results  for  BFSK  and  (120)  into  (114),  we  get  the 
total  probability  of  bit  error 


(2N-q')(2N-i-q\(  |  1 

*  l  2N  A  2N-\  As,,  +  2 J 


exp 


Further  simplification  results  by  using  the  substitutions 


which  implies 


and 


which  implies 


(122) 


(123) 


(124) 


(125) 


Then  using  (122)  through  (125)  in  (80)  and  (81) 


m 

f  Y  c+^cl 

o  hc^c 

lpexpi 

t  §c  J 

x  70 

L  V  4c  J 

(126) 


and 


( i ) 

f  yj+O 

-  \y Aj 

yexp. 

X  Iq 

[i  ^  J 

(127) 


b.  Special  Cases  when  Allowing  Two  Interference  Tones 

As  a  consequence  of  the  triple  nested  numerical  integration  required  by  (120) 
in  order  to  evaluate  (121),  obtaining  numerical  results  when  two  tones  can  jam  a  single 
hop  slot  is  computationally  intensive.  Based  on  the  results  when  thermal  noise  is 
neglected,  we  expect  the  performance  will  be  better  than  when  only  a  single  interfering 
tone  is  allowed. 

When  there  is  no  signal  fading 


yi  =  yb  =  C,b  and  =  0 

and  when  there  is  no  jammer  fading 

Tj  =  ^d  =  0 

In  the  case  of  Rayleigh  fading  of  the  signal 

fj,  =  4*  and  ^  =  0 


and 


yj  =  and  £/  =  0 


(128) 


(129) 


(130) 


(131) 


when  the  jammer  is  Rayleigh  faded.  Two  limiting  cases  allow  evaluation  of  the  most 
pessimistic  and  most  optimistic  performance.  Optimistically,  from  the  communicator's 
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point  of  view,  when  there  is  no  signal  fading  and  Rayleigh  fading  of  the  jammer,  the 
probability  of  bit  error  reduces  to 


Ph  = 


(2N-q\(2N-q-\\ 

vnr  )v~vrr)cxv 


-0 

+r?  V2^-?vxJ 

' 

2 

U^A2iv-i A  p 

[~y^  +  2J 

+ 


2 n)  (ifrf )  (i)  Jo*  JrAXy^)  x  ^[pNyj/q  ,  Jyb  +  2Nyj/q  +  2j 2Nybyj/q  cos8 


-[^]exp  y  (y b+  4Myj/q  +  2 J 2Nybyj/q  cos b) 


x  7o 


j2Nyj/q  J yb  +  INyjlq  +  2  ]2NybyJq  cos  8 


>dyjdb  (132) 


which  requires  two  numerical  integrations  rather  than  three. 

The  most  pessimistic  performance  result  from  Rayleigh  fading  of  the  signal 
and  no  fading  of  the  jammer.  In  this  case. 


n  (2N-q\(2N-q-\\ 

Pi  =  nFJn)v: TAxp 


A*] 

A* 

2 

UA2N-1  Axp 

[^  +  2j 

JrA^Y^}  *  { Q^Wyj/q  ,Jyh+  2Nyj/q  +  2  j2Nybyj/q  cos8 


exp 


y  (yb  +  4Nyj/q  +  2  j2Nybyj/q  cos  s) 


x  /o 


j2Ny.,/q  J yb  +  2Ny.,/q  +  2  JlNybyjlq  cos  8 


>dybdb  (133) 
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-(i)  exp  ■  (j)  (y*  +  ANyj/q  +  lj2Nyhyj/q  cos 


1 2Nyj/q  ,  Jyb  +  2Nyj/q  +  2  J2Nybyj/q  cos  5 


(134) 


where  a  single  finite  numerical  integration  over  the  phase  angle  is  required.  In  this  thesis 
only  ( 1 34)  is  evaluated.  Further  work  in  this  area  could  consider  the  effects  of  Ricean 
fading  in  the  both  tones  jammed  case  based  on  (120). 
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C.  MULTI-TONE  INTERFERENCE  NUMERICAL  PROCEDURE 


The  crux  of  the  numerical  computation  for  this  problem  is  the  accurate  and  efficient 
computation  of  Marcum’s  Q  function.  Originally  the  Q  function  described  a  shorthand 
notation  for  the  probability  integral  representing  the  output  of  a  correlation  detector 
containing  a  radar  target  in  the  presence  of  narrowband  Gaussian  noise  [Ref.  14]. 

The  Q  function,  described  by  (59),  has  no  closed  form  solution  expressible  in  a  finite 
number  of  terms.  Frequently  cited  equivalent  expressions  containing  an  infinite  series  of 
Bessel  functions,  while  valuable  as  analysis  tools,  do  not  ease  the  computational  burden. 
The  strategy  is  to  test  the  input  arguments  for  their  magnitude,  difference,  and  the 
presence  of  zeros.  Based  on  the  input  arguments  the  program  computes  the  value  of  the 
Q  function  in  an  appropriate  subroutine  and  returns  to  the  calling  program. 

First,  the  presence  of  a  zero  passed  as  an  input  argument  is  tested.  If  present  the 
following  simplification  results 

Q(a,  0)  =  1 .0  (135) 

and 

0(O,P)  =  exp{-(j)p2)  (136) 

The  second  consideration  is  the  magnitude  and  difference  of  the  input  arguments.  If  the 
product  of  the  arguments  is  greater  than  ten  and  their  difference  is  greater  than  five,  then 
the  Q  function  is  computed  using  an  asymptotic  polynomial.  If  the  magnitude  of  the 
product  is  greater  than  1000,  then  the  Q  function  is  computed  by  directly  integrating  (59) 
using  a  large  argument  approximation  for  the  Bessel  function.  For  the  semi-infinite 
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integration  required  here,  the  Romberg  technique  is  very  efficient.  The  Romberg 
technique  is  a  recursive  implementation  of  Richardson's  extrapolation  which  is  continued 
until  a  user  specified  tolerance  is  reached.  In  theory  and  in  practice  it  offers  much  faster 
convergence  and  smaller  errors  than  the  standard  trapezoidal  or  Simpson  technique. 

[Ref.  16] 

If  none  of  these  conditions  is  satisfied  the  program  defaults  to  a  computation  based 
on  a  numerical  inverse  Laplace  transform.  The  integrand  of  (59)  is  expressed  in  the 
Laplace  domain  as  [Ref.  10,  equation  6.643.4] 


TO=1-(^)xexpKfX'-2JTT 


(137) 


Then  the  integration  may  be  determined  from  the  inverse  transform  of 


Q(a,  P)  <=>  F(s)  =  \  x 


1  - 


2s 


Vt)* 


exp 


or 

2 


1  - 


25+1 


(138) 


This  is  useful  for  many  intermediate  values  of  the  input  arguments. 
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D.  MULTI  -TONE  INTERFERENCE  RESULTS 


First,  the  results  for  the  scenario  in  which  at  most  one  tone  is  jammed  per  hop  slot 
are  presented.  The  performance  for  a  wide  variety  of  channel  conditions  is  computed  for 
several  representative  signal  to  thermal  noise  ratios.  All  figures  are  for  one  hop  per  bit 
(L=  1 ).  For  each  channel  case  there  are  two  prime  questions  to  be  answered.  First,  "How 
does  performance  vary  for  fixed  values  of  qV  Second,  "Is  the  no  fading  analytic  and  no 
thermal  noise  worst  case  q  still  valid  for  fading  channels?"  To  answer  the  second 
question  the  proposed  worst  case  value  of  q ,  qQ  ,is  used  to  compute  performance  as  well 
as  some  values  of  q  offset  from  q0 .  To  show  the  trend,  the  results  of  q  +  4  and  q0  -4  are 
compared  with  qa  based  performance.  Table  2  provides  an  overview  of  the  channel 
conditions  considered. 
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TABLE  2.  MULTI-TONE  CHANNEL  CONDITIONS 


NAME 

Eb/N0  dB 

Rb  / 

Ricean  Fading 

Rj/ 

Ricean  Fading 

Figure  25 

13.35 

100 /Very  Weak 

100  /  Very  Weak 

Figure  26 

13.35 

100 /Very  Weak 

0.0  /  Rayleigh 

Figure  27 

13.35 

0.0  /  Rayleigh 

100  /  Very  Weak 

Figure  28 

13.35 

10  /  Moderate 

1 0  /  Moderate 

Figure  29 

13.35 

0.0  /  Rayleigh 

0.0/  Rayleigh 

Figure  30 

13.35 

1 0  /  Moderate 

1  /  Strong 

Figure  31 

13.35 

50  /  Weak 

50  /  Weak 

Figure  32 

13.35 

25 / Low 

25 / Low 

Figure  33 

20 

1 0  /  Moderate 

10  /  Moderate 

Figure  34 

20 

1  /  Strong 

1 0  /  Moderate 

Figure  35 

20 

1 0  /  Moderate 

1  /  Strong 

Initially,  we  consider  the  case  where  fading  has  very  little  influence  to  compare  with 
a  no  fading  analytic  result.  Figure  25  is  an  illustration  of  the  performance  obtained  with 
a  moderate  signal-to-thermal  noise  ratio  that  typically  provides  a  10'5  bit  error  ratio  in  the 
absence  of  fading  and  interference.  Thus,  Figure  25  provides  a  good  basis  for 
comparison  with  the  no  fading  case.  The  worst  case  performance  is  seen  as  the  envelope 
of  the  fixed  q  curves.  As  the  signal-to-jammer  power  increases,  the  performance  quickly 
approaches  2x10 5,  which  is  very  close  to  the  no  jamming  performance.  We  can  also  see 
that  the  degradation  inflicted  by  multi-tone  jamming  is  strongly  influenced  the  choice  of  a 
fixed  number  of  interference  tones. 
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Figure  25.  E,/N0=13.35dB,  R*  =100,  Rj=100 
To  bracket  the  range  of  expected  performance  the  case  of  extreme  of  Rayleigh 

jammer  fading  and  very  weak  signal  fading  is  depicted  in  Figure  26.  Here  the 

performance  is  slightly  improved  at  lower  signal-to-jammer  power  ratios  than  is 

observed  in  Figure  25.  Equivalently,  the  worst  case  envelope  is  slightly  inclined  more 

toward  the  y-axis  indicating  better  performance.  However,  at  low  signal-to-jammer 

power  levels,  the  performance  is  stiil  relatively  poor. 
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Figure  26.  E,/N0=13.35dB,  1^=100,  Rj=0.0 
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Figure  27.  E,/N0=13.35dB,  R*-  0,  Rj-100 
At  the  other  extreme,  a  pessimistic  performance  is  obtained  when  the  signal  suffers 

Rayleigh  fading,  but  the  jammer  has  very  weak  fading.  The  jammer  is  able  to  reduce 

performance  overall  and  maintain  a  high  bit  error  ratio  (BER)  for  higher  signal-to-jammer 

power  ratios.  Even  at  Eb/AO=30dB,  the  effects  of  fading  and  jamming  produce  a  BER 

near  4x1  O'2. 
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Figure  28.  E,/N0=13.35dB,  R^IO,  Rj=10 
In  Figure  28  we  display  the  performance  when  both  the  signal  and  the  jammer 

experience  moderate  Ricean  fading.  As  expected  the  performance  lies  between  the 

optimistic  case  of  Figure  26  and  the  pessimistic  case  of  Figure  27.  In  this  case  each  tone 

contains  the  same  direct-to-diffuse  power  ratio,  but  as  E JNj  grows  above  OdB,  the  direct 

signal  power  is  greater  than  the  direct  jammer  power  and  is  able  to  mitigate  the  influence 


of  multi-tone  jamming. 
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Figure  29.  E,/N0=13.35dB,  R^O,  Rj=0 

In  Figure  29,  the  signal  and  jammer  again  share  the  same  direct-to-diffuse  power 
ratio,  but  now  Rayleigh  fading  is  assumed.  The  performance  is  worse  overall  compared 
with  the  moderate  Ricean  case  and  approximately  3dB  worse  at  higher  signal-to-jammer 
power  levels.  Also  note  that  the  worst  case  occurs  at  q =1000.  Therefore,  the  greatest 
degradation  occurs  when  the  jammer  is  forced  into  the  multi-tone  equivalent  of 
broadband  jamming. 
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Figure  30.  E,/N0=i3.35dB,  R^IO,  Rj=l 

In  Figure  30,  the  channel  conditions  are  more  favorable  to  the  communicator  and 
represent  a  ten  fold  increase  in  the  Rb  over  the  Rj  from  the  performance  depicted  in 
Figure  28.  At  the  critical  value  of  OdB  for  the  signal-to-jammer  power  ratio  (EJNj)  the 
performance  in  the  favorable  channel  (  Figure  30)  is  3.3x1c4  while  in  the  moderate 
Ricean  channel  (Figure  28)  the  same  E JNj  produces  a  BER  of  2. lxl  O'3.  So  the  fact  that 
the  jammer  suffers  more  fading  over  the  Ricean  channel  produces  a  small  improvement 
for  the  communicator. 
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Figure  31.  E,/N0=l3.35dB,  R„=50,  Rj=50 
Figure  3 1  is  similar  to  Figure  25,  but  includes  less  fading.  The  performance  is 

essentially  the  same  as  the  very  weak  fading  performance,  but  with  some  slight 

improvement  for  EJNj  greater  than  OdB.  As  in  earlier  fading  channel  cases,  the 

performance  varies  greatly  for  fixed  values  of  q  when  the  jammer  has  a  power  advantage. 
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Figure  32.  E,/N0=13.35dB,  1^=25,  Rj=25 
Figure  32  depicts  another  drop  in  the  channel  direct-to-diffuse  ratio  shared  by  the 

signal  and  jammer.  At  this  level  the  difference  in  performance  below  the  very  weak  case 

where  the  direct-to-diffuse  ratio  was  100  is  noticeable.  When  the  jammer  enjoys  a  power 

advantage,  the  BER  is  still  unacceptably  high.  The  worst  case  envelope  curve  is  nearly 

linear  above  E,/No=0dB  and  nearly  flat  for  all  choices  of  q  above  OdB. 


77 


Pb 


Figure  33.  E,/No=20dB,  R^IO,  R,=10 

In  Figure  33  we  begin  to  look  at  the  performance  for  larger  signal-to-thermal  noise 
ratios.  An  increase  in  E,/N0  from  13.35dB  (Figure  28)  to  20dB  (Figure  33)  both  in 
moderate  Ricean  fading,  provides  an  overall  improvement  in  performance.  However,  the 
improvement  is  significant  when  EJNj  is  greater  than  -12dB  and  provides  nearly  a  factor 
of  ten  at  OdB. 
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Figure  34.  Eb/No=20dB,  Rt=10,  Rj=  10 

Figure  34  shows  how  rapidly  performance  can  degrade  even  for  a  strong  20dB 
signal-to-thermal  noise  power  ratio  when  the  signal  fading  progresses  from  moderate 
(Figure  33)  to  strong  Ricean  fading.  Similar  to  the  situation  observed  in  Figure  27  ,  the 
worst  case  performance  is  obtained  when  the  jammer  places  an  interfering  tone  in  every 
hop  slot.  The  jammer  therefore,  need  not  be  very  sophisticated  to  prevent  efficient 
communications.  In  Figure  35  the  roles  are  reversed.  With  the  same  signal-to-thermal 
noise  ratio,  the  signal  experiences  moderate  Ricean  fading  while  the  jammer  experiences 
strong  Ricean  fading. 


79 


Pb 


Figure  35.  E,/No=20dB,  R*=10,  Rj=l 

With  moderate  Ricean  signal  fading  a  particular  value  of  qa  is  required  to  cause  the 
worst  case  performance  demanding  greater  jammer  sophistication.  The  strongly  faded 
jamming  tone  is  much  less  effective  when  E JNj  is  greater  than  OdB.  Previously  we 
observed  an  approximate  tenfold  drop  in  BER  experienced  for  the  Eb/N0=13.35dB  case  as 
the  jammer  fading  grew  from  moderate  (Rj=  1 0)  to  strong  (R,=l )  at  Eb/7V,=0dB.  This 
trend  is  more  pronounced  at  the  Et/No=20dB  level  as  we  transition  from  Figure  34  to 
Figure  35  where  the  improvement  is  greater  than  a  factor  of  25. 
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The  influence  of  fading  on  the  selection  of  q  to  cause  worst  case  performance  in  the 
multi-tone  environment  remains  to  be  answered.  For  FFH/BFSK  in  the  absence  of  fading 
and  thermal  noise  we  expect  the  worst  case  to  occur  when  q  is  chosen  as  the  integer 
portion  of  the  jammer-to-signal  power  ratio.  A  comparison  of  the  following  worst  case 
performance  curves  will  show  them  to  be  the  envelope  of  the  previous  fixed  q  computed 
curves. 
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Figure  36.  Worst  Case  Multi-Tone  Jamming  for 
E,/N0  =  13.35dB.  R^O.  1^=100 

The  most  pessimistic  case  for  the  communicator  is  pictured  in  Figure  36  where  the 
signal  has  Rayleigh  fading  and  the  jammer  enjoys  essentially  no  fading.  W'e  see  that  in 
this  situation  the  jammer  does  not  need  to  be  very  sophisticated  to  hinder  communication. 
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The  results  displayed  in  Figure  36  and  in  Figure  37  indicate  there  is  essentially  no 
difference  in  the  q0  performance  and  the  q0  +4  and  q0  -4  cases. 
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Figure  37.  Worst  Case  Multi-Tone  Jamming 
Et/N0=13.35dB,  ^=0,^=0 

The  conclusion  to  draw  from  Figure  37  is  the  degradation  in  channel  conditions  for 
the  jammer  does  not  significantly  assist  the  communicator  when  both  are  acting  through 
Rayleigh  channels. 
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For  convenience  we  will  denote  the  q0+ 4  and  q0- 4  collectively  as  A q0.  It  is  not  until 
there  is  some  appreciable  direct  signal  power  reaching  the  receiver  that  the  q0  and  A q0 
begin  to  show  some  deviation. 

Pb 


Figure  38.  Worst  Case  Multi-Tone  Jamming 
Eb/No=20dB,  Rt,=10,  Rj=10, 

In  Figure  38  the  A q„  curves  dip  below  the  anticipated  worst  case  curve.  The  trend 
progresses  as  the  direct-to-diffuse  ratio  increases  in  Figure  39  and  Figure  40  respectively. 
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Figure  39.  Worst  Case  Multi-Tone  Jamming 
Eb/N0=13.35dB,  R*=25,  R,=25, 

Since  the  value  of  q  is  required  to  be  a  positive  integer  greater  than  or  equal  to  one, 
there  is  a  small  discontinuity  in  the  qa- 4  curve.  At  those  signal-to-jammer  power  values 
where  the  q0  value  would  be  less  than  5,  the  qQ- 4  defaults  to  <7=1  for  the  computation. 
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Figure  40.  Worst  Case  Multi-Tone  Jamming 
Eb/N0=13.35dB,  1^=50,  Rj=50, 


The  salient  observation  in  Figure  38,  Figure  39,  and  Figure  40  is  that  the  anticipated 
worst  case  performance  in  the  no  fading  and  no  thermal  noise  case  remains  the  worst  case 
performance  for  the  weakly  and  moderately  faded  channels  with  moderate  thermal  noise. 
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For  comparison  Figure  41,  Figure  42,  Figure  43  displays  the  composite  worst  case 


performance  for  various  signal-to-thermal  noise  ratios. 
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Eb/NJ  (dB) 


Figure  41 .  Worst  Case  Multi-Tone  Jamming  for  Several 
Values  of  E,/N0,  Rt=l,  Rj=l 

Figure  41  displays  the  performance  when  both  the  signal  and  the  jammer  suffer 
strong  fading.  The  diffuse  nature  of  the  channel  keeps  the  probability  of  bit  error 
relatively  high  despite  high  signal-to-jammer  power  ratios  and  serves  more  to  the 
detriment  to  the  communicator  than  the  jammer. 
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When  the  channel  has  moderate  vice  strong  Ricean  fading  the  performance  is  that 
shown  in  Figure  42.  Above  -6dB  EJNj  there  is  a  significant  improvement  in 
performance  over  that  in  Figure  41.  The  noticeable  change  in  slope  for  the  20dB  curve  is 

Pb 


Figure  42.  Worst  Case  Multi-Tone  Jamming  for  Several 
Values  of  E,/N0 ,  R,,=  1 0,  R;=  1 0 

a  consequence  of  the  signal-to-jammer  phase  term  in  the  bit  error  analysis.  When  the 
signal  amplitude  is  consistently  greater  than  the  jammer's  amplitude,  the  phase  difference 
has  little  effect  and  the  rate  of  improving  performance  increases.  This  slope  change  is 
also  noticeable  in  Figure  43  for  the  20  dB  curve. 
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Figure  43.  Worst  Case  Multi-Tone  Jamming  for  Several 
Values  of  Et/N,,,  1^=100,  Rj=l 

Figure  43  is  a  composite  of  performance  under  channel  conditions  favorable  to  the 
communicator.  However,  the  10:1  ratio  in  the  signal's  direct-to-diffuse  over  the  jammer's 
direct-to-diffuse  ratio  provides  little  relief  to  the  communicator. 

Figure  44  displays  the  performance  when  the  jammer’s  intelligence  is  relaxed  to 
allow  two  interference  tones  in  a  hop  slot  when  neither  the  signal  nor  the  jammer 
experience  fading. 
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Figure  44.  Worst  Case  Multi-Tone  Interference  when 
Allowing  Two  Jamming  Tones  per  Hop  Slot 

Figure  44  exhibits  the  same  linear  improvement  in  performance  as  EJNj  decreases 
for  several  values  of  E,/N0 .  There  is  a  marked  change  in  slope  at  E b/N,  =0dB.  Overall 
performance  is  better  thar.  the  perffomance  when  at  most  one  tone  per  hop  slot  is  jammed 
in  a  near  Gaussian  channel  as  depicted  in  Figure  40. 
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IV.  CONCLUSIONS 


A.  FFH/MFSK  PARTIAL-BAND  NOISE  JAMMING 

A  performance  analysis  has  been  completed  for  the  conventional  fast 
frequency-hopped  M-ary  orthogonal  frequency-shift  keyed  receiver  employing  diversity 
with  noncoherent  reception  of  signals  transmitted  over  Ricean  fading  channels.  We  can 
draw  several  conclusions  based  on  the  results. 

Diversity  alone  is  insufficient  to  overcome  the  effects  of  worst  case  partial  band 
jamming  acting  in  Ricean  fading  channels.  The  intelligent  jammer  can  optimize  the 
fraction  of  bandwidth  jammed  to  force  a  worst  case  bit  error  ratio  (BER).  The 
noncoherent  combining  losses  experienced  by  noncoherent  systems  is  aggravated  by  the 
low  signal-to-total  noise  ratio.  These  two  factors  combine  to  mitigate  the  advantages  of 
increased  diversity  in  a  constant  energy  per  hop  system. 

Increasing  the  modulation  order  provides  only  modest  improvements  in  performance 
which  decrease  as  diversity  increases.  For  most  channel  conditions,  the  performance  for 
A/=4  and  A/=8  is  essentially  the  same  for  all  L  greater  than  four.  The  only  remaining 
advantage  for  increased  modulation  order  is  to  offset  the  reduction  in  data  rate  that  results 
from  increasing  diversity,  L,  in  a  constant  energy  per  hop  system. 

When  the  signal  experiences  Rayleigh  or  strong  Ricean  fading,  the  jammer  need  not 
be  very  sophisticated  to  cause  marked  impairment  to  the  communicator.  Additionally, 
when  comparing  the  partial-band  jamming  strategy  for  A/=2  and  L=  1  in  which  the  signal 
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suffers  moderate  to  strong  fading  to  multi-tone  jamming  under  the  same  conditions,  we 
see  that  the  partial-band  jamming  strategy  causes  a  higher  BER  while  requiring  less 
precise  know  ledge  of  the  receiver  by  the  jammer. 

B.  FFH/BFSK  MULTI-TONE  JAMMING 

In  this  thesis,  the  case  of  l.~\  is  considered  for  the  fast  frequency -hopped  receiver 
described  above  The  influence  of  Ricean  fading  of  both  the  signal  and  the  lammcr  is 
included  in  the  analysis  as  is  the  effect  of  thermal  noise 

first,  the  counter-intuitive  result  is  obtained  that  Ricean  fading  of  the  lammer 
provides  little,  if  any.  relief  to  the  communicator,  while,  as  expected,  numerate  to  strong 
fading  of  the  signal  tone  marked  I  \  reduces  performance  Worst  case  multi-tone  lamming 
combined  with  signal  fading  produces  unacceptably  high  probabilities  of  bit  error,  i  e  . 
greater  than  10  Above  this  level  it  is  unlikely  that  forward  error  correction  coding  can 
be  successfully  implemented 

The  intelligent  tammer  can  select  the  number  of  tones  to  jam  which  pnnluces  a  worst 
case  performance  The  value  of  if  chosen  in  the  fading  channel  with  thermal  noise  to 
cause  worst  case  performance  is  the  same  value  of  if  armed  at  analy  tically  for  the  noise 
free  case  w  ithout  fading  For  all  combinations  of  channel  conditions,  not  until  the 
signal-to-jammcr  power  ratio  exceeds  lOdB  is  the  performance  essentially  thermal  noise 
limited  When  the  jammer  has  a  3dB  or  better  power  advantage,  the  performance,  in  cv  en 
the  most  optimistic  channel  conditions,  is  very  poor  and  is  dominated  by  the  tammer 
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